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Abstract 

Sulfuric acid (H2SO4), formed in commercial aircraft operations via fuel— S — > SO2 — > SO3 — > 
H2SO4, plays an important role in the formation of contrails. It is believed that the first step 
occurs inside the combustor, the second step in the engine exit nozzle, and the third step in the 
exhaust plume. Thus, measurements of the sulfur oxidation rates are critical to the understanding 
of contrail formation. Field measurements of contrails formed behind commercial aircraft 
indicate that significantly greater conversion of fuel-bound sulfur to sulfate aerosol occurs than 
can be explained by our current knowledge of contrail physics and chemistry. The conversion of 
sulfur from S(IV) to S(VI) oxidation state, required for sulfate aerosol formation, is 
thermodynamically favored for the conditions that exist within jet engines but is kinetically 
disfavored. The principal reaction pathway is O+SO2+M — > SO3+M. The rates of this reaction 
have never been measured in the temperature and pressure regimes available to aircraft 
operation. 

In the first year (FY 02 ) of this project, we performed a series of experiments to elucidate the rate 
information for the O+SO2+M — » SO3+M reaction. The work performed is described below, 
following the proposed work plan. Because we used the H2/O2 system for an O-atom source and 
rate coefficients were obtained via computer simulation, construction of a reaction mechanism 
and either recalculation or estimation of thermodynamic properties of H x SO y species are 
described first. 


Reaction Mechanism of the SO2/H2/O2 System 

For the H2/O2 system, we adopted a reaction mechanism of 20 elementary reactions that has been 
tested extensively at various conditions. 1 ” 3 Then, additional reactions for SO2 oxidation by H x O y 
species in the SO2/H2/O2 system and their rate coefficients were critically reviewed and a 
reaction mechanism of 19 reactions for this system was composed. A comprehensive but 
tnincated reaction mechanism for the SO2/H2/O2 system is given in Table I of Appendix A. A 
brief review of important rate coefficients is given in the Discussion section. 

Thermodynamic data of relevant species of H2/O2/SO2 System 

In computer modeling using a reaction mechanism for practical combustion systems or for 
extracting rate coefficient information, reliable thermodynamic properties (A fF/°( 298. 15 K), 
H(T), S(T), G(T), C P (T )) of the reacting species are necessary. Recently the value of Z/WisfH- 
OH)=l 1 8.81 ± 0.07 kca\ mol” 1 was measured by Ruscic et al. 4 From this new bond dissociation 
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enthalpy of H 2 O, Af// J (298.15 K) of OH=8.91 ± 0.07 A:cal mol -1 has been reestablished and the 
differences from the widely used value of Gurvich et al . 5 or JANAF 6 have been discussed in 
great detail . 4 Previously we had used 9.4 1 ± 0.05 /teal mol ' 1 for A f/7°(298. 1 5 K) of OH, based 
upon the recommendation of Gurvich et al. The difference of 0.49 kcal mol -1 from the new value 
of Ruscic et al. is quite large. Considering the uncertainty in obtaining D 0 (O-H) by Gurvich et al. 
from spectroscopic measurements and an unambiguous direct measurement of Z)°(H-OH) by a 
photoionization experiment by Ruscic et al., we adopted the new AfH°(298.15 K) of 
OH=8.91 ± 0.07 /teal mol -1 of Ruscic et al. and readjusted all thermodynamic properties of OH 
accordingly. The thermodynamic properties of other species in the H 2 /O 2 system are unchanged. 

While thermodynamic properties of species in the H 2 /O 2 system are more or less well 
established, those of H x SO y species (HS, HSO, HOSO, HSO 2 , and HOSO 2 ) are virtually 
unknown. Using available experimental data and theoretical calculations, thermodynamic data of 
these species are estimated. Details are described in Appendix B. 

Selection of Experimental Conditions via Sensitivity Analysis 

In determining rate coefficients of the reaction 

O+SO 2 +M SO 3 +M (R21) 

an O-atom source must be found. In this work O-atoms were produced from H 2 /O 2 reactions. I.e., 
the H 2 /O 2 reaction system was perturbed with a small amount of SO 2 . The major reactions that 
consume and produce O-atoms in the H 2 /O 2 system are, respectively 


0+H 2 OH+H (R2) 

and 

H+0 2 OH+O (Rl) 


Hence, in the SO 2 /H 2 /O 2 system, R21 would compete with R2. If we use rich H 2 /O 2 mixtures, O- 
atoms would be consumed mainly by R2. This restricts us to consider lean H 2 /O 2 mixtures only. 
Figure 1 shows R=|log([/] P k/[H 2 ]o)| as a function of equivalence ratio <|) of H 2 /O 2 mixtures, where 
/= O, OH, and H. Small R values indicate production of a large amount of species /. Obviously, 
lean mixtures are more suitable than rich mixtures in this regard. However, at given jet engine 
exhaust nozzle temperatures (900 to 1200 K), the following reaction will compete with Rl for H- 
atoms 


H+O 2 +M HO 2 +M (R9) 

Thus, it is also desirable to select reaction mixtures that produce a less amount of H-atoms. 

Sensitivity analyses for an elapsed time to reach half of the absorption maximum ( T 50 ) were 
performed using the following baseline H 2 /O 2 mixture (without SO 2 ) with various levels of 
sensitization by SO 2 at T=1000 K, p =1 5.0 pmol cm -3 ; 0.5% H 2 / 10 . 0 % O 2 (<|)=0.025), 0.15%- 
3.0% SO 2 / 0.5% H 2 / 10 . 0 % O 2 (0 =0.025)/ Ar. In the baseline mixture, R9, Rl, R2, and R3 show 
sensitivities for T 50 in that order. Figure 2(a) shows a sensitivity spectra for 0.5% H 2 / 10 . 0 % 
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0; : /89.5% Ar mixture. Adding 0.15%, 0.25%, and 0.5% SO 2 to a 0.5% H2/10.0% Ch/Ar mixture 
produced only a trace level of sensitivity for R2 1 . The sensitivity of R2 1 was found to be 
dependent upon the SO 2 /H 2 ratio in the mixture. Furthermore, as the SO 2 /H 2 ratio is increased, 
the sensitivity of R21 and R9 is notably increased together with the relative importance of R21 
over R9. Figure 2(b) shows a sensitivity spectra for the 2.0% SO 2 / 0.5% Fh/10.0% 02/88.5% Ar 
mixture. Based upon this sensitivity analyses we selected the following three mixtures: for the 
baseline experiments, 0.5% H2/10.0% 02/89.5% Ar; for sensitization experiments, 2.0% 
SO 2 /0.5% H 2 /10.0% 0 2 /87.5% Ar and 3.0% SO 2 /0.5% H 2 /10.0% 0 2 /86.5% Ar. 

Experimental 

The shock tube-laser absorption technique is described in detail elsewhere. 3 Briefly, the test 
section of the shock tube was routinely pumped below 3 p Torr by a Varian V60 Turbopump. 

The combined leak and out-gassing rate of the test section was 5 p Torr min -1 . All shocks were 
initiated within 1 minute of admitting the test gas mixture. Shock velocities were measured using 
4 flush-mounted 1 13A21 PCB Piezotronics pressure transducers coupled to Phillips PM6666 
programmable counter-timers. The last pressure transducer was axially coincident with the center 
of the observation windows, 12.7 mm distant from the end wall. Shock speeds were fitted to a 
second-order polynomial in distance. The extrapolated speed to the end wall was used together 
with NASA thermodynamic data 7 in solving the standard relationships for the reflected shock 
properties, assuming full vibrational relaxation and no chemical reaction at the shock front. The 
computed reflected shock properties were corrected for boundary layer interaction effects using 
the similar method of Michael and Sutherland. 3,8,9 

The combined OH and SO 2 absorption profiles were measured at 310.032 nm (in air), 
corresponding to the PI (5) line of the (0,0) band of the OH A 2 Z + X 2 n transition — at this 
wavelength, SO 2 also shows some absorption. The UV light was produced by a CW argon ion 
laser (Coherent Innova 200) pumped ring-dye laser (Coherent 899 - 21) running Kiton Red 620 
dye. Intracavity frequency doubling was achieved with an angle tuned LiIC >3 crystal. For the 
detection of OH radicals and SO 2 molecules, a triple beam scheme was employed; probe beam, 
intensity reference beam and a wavelength reference beam (through a burner-stabilized 
methane/air flame). The detector was a combination of interference filter (3 1 0 nm with 1 0 nm 
fwhm) — PMT(Thom EMI9924QB) — high speed buffer/amplifier. The S/N of the "probe beam- 
intensity reference beam" was usually better than 100. The overall electronic time constant 
determined for the entire PMT-buffer-cable system was less than 0.5 p s. 

H 2 / 02 /Ar and S 02 /H 2 / 02 /Ar test gas mixtures were prepared manometrically and allowed to 
stand for 48 hr before use. Gases were used without further purification. The stated purities were: 
H 2 , 99.9995% (MG Industries, Scientific Grade); O 2 , 99.998% (MG Industries, Scientific 
Grade); SO 2 (anhydrous), 99.98 wt.% (MG Industries, Scientific Grade); Ar, 99.9999% (MG 
Industries, Sputtering Grade). 
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Results 


Because 310 nm UV light is also absorbed by SO2 molecules, it is necessary to determine 
absorption coefficients of SO2 ( e e (S02). The values of 8 e (S02> determined are shown in Figure 
3 . A least squares fit to the data gives an expression 

8 e (S0 2 )=2.16x 10 5 exp(-765.0 K/T) cm 2 mol -1 

Evaluation of a characteristic time ( x 50) in the baseline experiments was done in the usual way. 
However, in the sensitization experiments, because there have always been long delays for 
sensible absorption signal rise (OH production), initial absorption by SO2 immediately after the 
reflected shock was subtracted from the total absorption signal and x 50 values were evaluated as 
in the baseline experiments. Figure 4 shows X50 values measured for the baseline mixture, 0 . 5 % 
H2/10.0% 02 / 89 . 5 % Ar, and for the sensitization mixtures, 2 . 0 , 3 . 0 % SO2/0.5% H2/10.0% 
O2/86.5, 87 . 5 % Ar. Throughout all of the experiments, densities were kept nearly the same 
between the baseline and sensitization experiments. 

As seen in the figure, the X 50 values of sensitization mixtures are longer at low temperatures and 
increase with SO2/H2 ratio. At temperatures above 1200 K, the differences disappear between the 
two levels of sensitization. 


Discussion 


H 2 /O 2 System 

Because we generated O-atoms from the H 2 /O 2 reaction system, a complete understanding of the 
H 2 /O 2 reaction chemistry is a prerequisite for rate coefficient determination of the O+SO 2 +M — > 
SO 3 +M reaction (R21). There are three key reactions in H 2 /O 2 chemistry: chain branching H+O 2 
— » OH+O (Rl), chain propagating O+H 2 —> OH+H (R2) and chain terminating H+O 2 +M — > 
HO 2 +M (R9). A comprehensive review of these three reactions is given here. 

(a) Chain Branching H+O 2 — » OH+O (Rl) 

Figure A(l) in Appendix A shows the rate coefficient expressions reported recently. The current 
widely used rate coefficient expressions are those of Du and Hessler , 10 Baulch et al . 11 and GRI . 12 
Du and Hessler obtained their kj expression, £j=9.76x 10 13 exp(-7474 K/T) cm 3 mol -1 s _1 (960 to 
5300 K) by combining their rate coefficient measurements at very high temperatures (3450 to 
5300 K) with the results of Masten et al . 13 (1450 to 3370 K), Shin and Michael 14 (1 103 to 2055 
K) and Pirraglia et al . 15 (962 to 1 705 K). The compilation by Baulch et al. gives an expression 
identical to that of Du and Hessler. Pirraglia et al. combined their relatively low temperature 
results with the high temperature data of Frank and Just 16 (1693 to 2577 K) and reported 
/t,=1.9xl0 14 exp(-8272 K/T) cm 3 mol ” 1 s _1 (962 to 2577 K). The GRI expression, £i=2.65xl0 16 
T 0 6707 exp(-8575 K/T) cm 3 mol ” 1 s " 1 (300 to 3000 K) is a fit to the combined high temperature 
results of Yu et al . 17 (which, in turn, the results of Yuan et al . 18 and of Masten et al.) and the 
somewhat low rate coefficients of Pirraglia et al. (especially low temperature data) with the room 
temperature data compiled by Demore et al . 19 Previously, we (GRC ) 1 reported &i=7.13x 10 13 
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exp(-6957 K/T) cm 3 mol 's 1 (1050 to 2500 K) based upon OH radical growth rate that 
eliminates any H-atom producing contaminant effect. 

Since the Af//°(298.15 K) value of OH has been remeasured (see above), virtually all k\ values 
must be reexamined. Hence, with the new thermodata of OH we reevaluated our existing data. 
The rederived expression is £i=3.74x 10 15 T ~ 0467 exp(-7758 K/T) cm 3 moT 1 s _1 . Extrapolation 
of this expression to room temperature gives about a factor of 2.5 higher than the value 
recommended by Demore et al. As seen, this expression has a negative temperature dependence 
of the A-factor. The temperature dependence of the A-factor has been given by several authors: 
-0.907 by Schott , 20 -0.927 by Yuan et al., 18 -0.7 by Masten et al . 13 and -0.816 and +0.500 by 
theoretical calculations of Miller 21 and Miller and Garrett , 22 respectively. 

Due to the difficulty in performing experiments at relatively low temperatures, no group has 
extended their temperature range wide enough to obtain a consistent set of experimental data. 
Instead, the results of Pirraglia et al. (960 to 1200 K) have been extensively used for the fit of the 
k\ expression by various authors. Recently we measured OH absorption profiles in various H 2 /O 2 
systems at temperatures down to 950 K and up to 3500 K. We are currently analyzing these 
results. A combination of these new results with the previous ones reported in this laboratory 1 
will form a consistent set of data over the widest temperature range (950 to 3500 K) ever 
employed by a single laboratory. Based upon this data set, a new rate coefficient expression will 
be produced and the validity of its temperature dependence of the A-factor and the extrapolated 
value to room temperature will be tested. 

(b) Chain Propagating O+H 2 — > OH+H (R2) 

SO 2 molecules in the SO 2 /H 2 /O 2 system are consumed mainly by O+SO 2 +M —> SO 3 +M (R21). 
The chain propagating reaction R2 inherently competes with R21 for O-atoms. Figure A(2) in 
Appendix A shows representative k 2 expressions published. Previously we presented 
kr=\.S7x 10 14 exp(-6854 K/T) cm 3 moT 1 s~’ (1420 to 2430 K). Reevaluation of OH absorption 
data obtained for k 2 evaluation with the new OH thermodata yielded k 2 -2.0A'/, 10 14 exp(-6972 
K/T) cm 3 moT 1 s" 1 . The differences between the old and new expressions are less than 3.5% in 
the given temperature range, the new expression being higher with temperature. One notable 
expression for k 2 is that of Sutherland et al . 23 Sutherland et al. reported a non-Arrhenius type rate 
coefficient expression, k 2 =5.06x 10 4 T 2 67 exp(-3480 K/T) cm 3 moT 1 s~ ! by combining their 
results of flash photolysis-shock tube work (880 to 2495 K) and the flash photolysis-resonance 
fluorescence technique (504 to 923 K) with the low temperature data of Presser and Gordon (297 
to 471 K ). 24 Our k 2 expression represents the high temperature data (880 to 2495 K) of 
Sutherland et al. quite well. However, due to the nature of a non-Arrhenius type curve fit, the 
extended temperature fit of Sutherland et al. over-predicts their results at temperatures around 
950 K and above 2500 K. Therefore, in our temperature range (950 to 1100 K), extrapolation of 
our expression should be more appropriate. 

(c) Chain Terminating H+O 2 +M — > HO 2 +M (R9) 

The chain terminating R9 reaction plays a governing role for ignition and flame characteristics at 
high pressure and especially at low temperature conditions (see sensitivity spectra in Figure 2). 
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For this reason, many authors 1, 25-42 investigated this reaction to determine rate coefficients. 
Experimental data and k 9 expressions are plotted in Figure A(3). At room temperature a 
consensus value has been established . 33 However, in the temperature range of combustion 
(T>750 K), the scatter of rate coefficient data is very large. For example, at T=~1000 K, there is 
a factor of 3 scatter in the £9 values. We performed a series of experiments to determine the k 9 
values in the temperature range of 950 to 1 150 K using a shock tube — laser OH absorption 
spectroscopy technique in H 2 /O 2 system. Experimental data were fit to £ 9 = 3.75 x 10 17 T 0 7 cm 6 
mol -2 s -1 . Our values positioned somewhat in the middle of the data scatter at T=~1000 K. 
Extrapolation of our kg expression also well represents the data of Hsu et al . 32 at temperatures 
between 296 K and 640 K and the room temperature result of Carleton et al . 33 

SO 2 /H 2 /O 2 System 

As mentioned before, sensitization was used with 2 % and 3% SO 2 . Sensitization resulted in a 
remarkable increase of t 50 (see Figure 4). The effect of sensitization level on x 50 becomes more 
noticeable as temperature decreases, but it diminishes as temperature increases. Nevertheless, 
this increase of x 50 is due to depletion of O-atoms by R21, so the H-atom production rate by R2 
is reduced. This, in turn, reduces the branching rate of Rl. 

The rate coefficients of R21 are derived by computer simulation using the Table I mechanism in 
Appendix A. In computer simulation, at each experimental condition, only the starting £21 value 
was varied so as to minimize the difference between the computed and experimental T 50 . Figure 
5 shows the provisional £21 values obtained in this work. R21 is a pressure (density) dependent 
reaction. For this type of reaction, it is usual to use two rate coefficient expressions at the 
extreme conditions, a low pressure limit fei.o and a high pressure limit £ 21,00 The rate 
coefficients in the fall-off region are then obtained by interpolation using these two extreme rate 
coefficients. Obviously R21 is in the falloff region in our experimental conditions. Therefore we 
computed the starting £21 value using the £ 21 , 0 , £ 21 , 00 , and a simple Lindemann-Hinshelwood 
interpolation formula 


£2 1 ^£2 1 , 00 — (£2 1 ,o /£2 1 , 00 )/( 1 +£2 1 ,q /^2 1 , 00 ) = x/( 1 +x) 


where x=£ 2 i,o/£ 2 i, co 


For £21,00 and £ 21 , 0 . we used the expressions given by Mueller et al ., 43 £ 21 , 00 “9.2 x 10 10 exp (-1200 
K/T) cm 3 mol -1 s ’ 1 and by Troe , 44 £ 2 i,o = 4 . 0 x 10 28 T -40 exp(-2650 K/T) cm 6 mol -2 s -1 , 
respectively. The £21.00 of Mueller et al. is the only high pressure limit expression ever reported. 
Several measurements of £ 21,0 have been reported either at high temperatures 45,46 (1435 to 2500 
K) or at low temperatures (<430 K ). 19,47-50 In the temperature range of our interest (900 to 1 100 
K), no experimental measurements are available. The £ 21,0 expression shown above was made by 
Troe by fitting the low and high temperature data. Note that at T<500 K, the rate coefficient 
value increases with temperature (positive activation energy), while at T>500 K, it decreases 
with temperature (negative activation energy). Our provisional £21 values obtained in the falloff 
region are about a factor of 2 to 3 higher than the theoretical calculations using the equation 
shown above. 
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Conclusions 


Existing data for H 2 /O 2 and SO 2 /H 2 /O 2 chemistry were critically reviewed and a reaction 
mechanism was assembled (Table I in Appendix A). Optimum H2/02/Ar and S02/H2/02/Ar 
mixture compositions and experimental conditions were established via a series of sensitivity 
analyses. Shock tube experiments were performed using OH laser absorption spectroscopy for 
the H 2 / 02 /Ar mixtures to examine baseline behavior and for the S 0 2 /H 2 / 02 /Ar mixtures to 
determine &21 values for the conditions selected. The experimental results were analyzed via 
computer simulations and the provisional rate coefficients of R21 were obtained. The provisional 
rate coefficients of R21 (£21 values) obtained in our experimental conditions are about a factor of 
2 to 3 higher than the theoretical calculations using the expressions of 1,0, £21,°° , and the 
parameters of the interpolation formula in Table I of Appendix A. Additionally, temperature 
dependent absorption coefficients of SO 2 at OH absorption wavelength were measured. 

Measured or calculated thermochemical and structural data for HS, HSO, HOSO, HSO 2 , and 
HOSO 2 species were critically reviewed and thermodynamic data for these species were obtained 
(Appendix B). 


Future Work 

In the second year of this project, besides what is contained in the project plan, the collision 
efficiency effect of SO 2 the on pressure dependent reactions, R9 and R2 1 , and the interpolation 
formula used in determination of ki\ will be examined further. In pressure dependent reactions, 
reaction rates increase with pressure. The increase of collision efficiencies of collision partners is 
equivalent to an increase of density (pressure). If one assigns a high collision efficiency of SO 2 
to R9 and R21, their reaction rates and x 50 values will increase and result in a lowering of ki\ 
values. (Currently we used a collision efficiency of 1 for SO 2 for R9 and R21 .) Experiments at 
various densities will also be performed to optimize the falloff parameters. Data analysis will be 
completed for the experiments performed to determine the rate coefficients of the chain 
branching reaction Rl. 
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|Log([i l pk /[H 2 lo)| vs. <(>; 1=0, OH and H 
T - 1000 K, p « 15.0 nmol cm 3 



Log(<|») 


Figure 1. — Absolute ratios of the peak concentrations of O- and H-atoms and OH radicals to the initial concentration 
of H 2 in the H 2 /0 2 /Ar mixtures as a function of the mixture equivalence ratio, <|) , evaluated at T=1000 K and 
p =15.0JImol cm 3 . Large values indicate the suppression of peak concentrations. 


0.5% H 2 , 10.0% 0 2 , 89.5% Ar 
T=1000 K, p-15.0 fimol cm 4 

h 2 +o 2 +m==ho 2 +m 
OH+H 2 =H 2 OHH 

g+h 2 =oh+h 
h+o 2 =oh+o 

-1 - 0.5 0 0.5 1 1.5 

Sensitivity 

Figure 2(a). — Sensitivity spectra for X 50 of the baseline mixture, 0.5% H 2 /10.0% 0 2 /89.5% Ar. 

T=1000 K and p =15.0jLl mol cm 3 . 
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Ln[e,(S0 2 )/cm 2 mol 1 ] 


2.0% SOj, 0.5% H 2 , 10.0% Oj, 87.5% Ar 
T=1000 K, p=15.0 |xmol cm -1 



Sensitivity 

Figure 2(b). — Sensitivity spectra for X 50 of the 2.0% SO 2 /0.5% H 2 /10.0% 0 2 /87.5% Ar. 
T=1000 K and p =15.0/1 mol cm -3 . 


Absorption Coefficients of S0 2 at 310 nm 
£ e (S0 2 )=2.16 x 10 5 exp(-765.0 K JT) cm 2 mol 1 


13 
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12 

11.5 
11 

10.5 - 

7.0 8.0 9.0 10.0 11.0 

10000 KIT 

Figure 3. — Base ‘e’ absorption coefficients of S0 2 at 310 nm. The least squares fit 
expression is £ C (S0 2 )=2.16X 10 5 exp(-765.0/T) cm 2 moT 1 . 
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Log(£ 21) ^ Logdso/jis) 


S0 2 /H 2 /0 2 /Ar 
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9.0 9.5 10.0 10.5 

10000 KIT 

Figure 4. — Measured T 50 values of the baseline mixture , 0.5% H 2 /10.0% 0 2 /89.5% 
Ar and two sensitization mixtures, 2.0% SO 2 /0.5% H 2 /10.0% 0 2 /87.5% 

Ar and 3.0% SO 2 /0.5% H 2 /10.0% 0 2 /86.5% Ar. 
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O SOj/H 2 /0,/Ar=2.0/0.5/10.0/87.5, p=16.1 pmol cm' 3 
□ S0 2 /H 2 /0,/Ar=3. 0/0. 5/1 0.0/86. 5, p=15.9 pmol cm' 3 
A S0 2 /H 2 /0 2 /Ar=0.0/0.5/1 0.0/89.5, p=15.4 pmol cm' 3 


O + SO, (+M) = SO3 (+M) 



Figure 5. — The falloff k 2 \ values of the reaction R21, 0+S0 2 (+M) — > S0 3 (+M). 
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Appendix A 


Table I. — Reaction Mechanism 0 



reaction 

A 

n 

© 

Note 

1 

h+o 2 =oh+o 

7.13E+13 

0.0 

6957 


2 

o+h 2 =oh+h 

1.87E+14 

0.0 

6854 


3 

oh+h 2 =h+h 2 o 

2.16E+08 

1.51 

1726 


4 

o+h 2 o=oh+oh 

4.51E+04 

2.70 

732 


5 

o+o+m=o 2 +m 

1.0E+17 

-1.0 

0 



Ar, 1 .0/ H 2 ,2.9/ 0 2 , 1 .2/ H 2 0, 1 8.5/ 



e 

6 

h+h+m=h 2 +m 

6.4E+17 

-1.0 

0 



Ar,1.0/ H 2 ,4.0/ H 2 O,12.0/ H,26.0/ 



e 

7 

H+0+M=0H+M 

6.2E+16 

-0.6 

0 



Ar,1.0/ H 2 O,5.0/ 




e 

8 

h+oh+m=h 2 o+m 

8.4E+21 

-2.0 

0 



Ar,1.0/ H 2 ,2.5/ H 2 0, 16.25/ 



e 

9 

h+o 2 +m=ho 2 +m 

7.55E+17 

-0.8 

0 



Ar,1.0/ 0 2 ,l-33/ H 2 ,3.33/ H 2 0,21.3/ 



e 

10 

ho 2 +h=oh+oh 

1.5E+14 

0.0 

505 


11 

ho 2 +h=h 2 +o 2 

2.5E+13 

0.0 

350 


12 

ho 2 +h=h 2 o+o 

5.0E+12 

0.0 

710 


13 

ho 2 +o=o 2 +oh 

2.0E+13 

0.0 

0 


14 

ho 2 +oh=h 2 o+o 2 

2.0E+13 

0.0 

0 


15 

H02^"H02 = 02^ - H 2 0 2 

1.3E+11 

0.0 

-820 

f 


H02+H02“02+H 2 02 

4.2E+14 

0.0 

6039 

f 

16 

h 2 o 2 +m=oh+oh+m 

1.2E+17 

0.0 

22900 



Ar,0.67/ O 2 ,0.78 / H 2 O,6.0/ 



e 

17 

h 2 o 2 +h=ho 2 +h 2 

1.7E+12 

0.0 

1900 


18 

h 2 o 2 +h=h 2 o+oh 

1.0E+13 

0.0 

1805 


19 

h 2 o 2 +o=ho 2 +oh 

2.8E+13 

0.0 

3225 


20 

H 2 0 2 +0H=H0 2 +H20 

7.0E+12 

0.0 

720 


21 

S0 2 +0+M— SO3+M 

9.2E+10 

0.0 

1198 

b 



4.0E+28 

-4.0 

2642 

c 


1.000 

1.000 

0.0 

0 

d 


Ar,1.0/ 0 2 ,1.33/ SO 2 ,1.0/ 




e 

22 

so 2 +oh+m=hoso 2 +m 

7.2E+12 

0.0 

360 

b 



4.5E+25 

-3.3 

360 

c 


0.640 

OO 

300.0 

0 

d 


Ar,1.0/ 0 2 ,1.5/ SOjMAf 




e 

23 

hoso+m=so 2 +h+m 

1.7E+10 

0.8 

23602 

b 



1.6E+31 

-4.53 

24759 

c 


0.450 

0.0 

OO 

OO 

d 


Ar,0.75/ H 2 ,2.5/ H 2 O,12.0/ 



e 

24 

so 2 +oh=so 3 +h 

4.9E+02 

2.69 

11977 
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Table I. — Reaction Mechanism 0 (continued) 



reaction 

A 

n 

0 

Note 

25 

so+o+m=so 2 +m 

3.2E+13 

0.0 

0 

b 



2.9E+24 

-2.9 

0 

c 


0.550 

0.0 

00 

OO 

d 


Ar,1.00/ 02,1.33/ H 2 O,10.0/ 



e 

26 

SO3+O— so 2 +o 2 

4.11E+11 

0.0 

3070 


27 

so 3 +so=so 2 +so 2 

1.0E+12 

0.0 

2013 


28 

so 3 +h 2 o=h 2 so 4 

7.23E+08 

0.0 

0 


29 

SO+OH+M=HOSO+M 

8.0E+21 

-2.16 

418 



Ar,0.75/ H 2 ,2.5/ H 2 O,12.0/ 




e 

30 

so+oh=so 2 +h 

5.2E+13 

0.0 

0 


31 

so+o 2 =so 2 +o 

6.2E+03 

2.42 

1535 


32 

hoso+h=so+h 2 o 

6.3E-10 

6.29 

-956 


33 

hoso+oh=so 2 +h 2 o 

1.0E+12 

0.0 

0 


34 

hoso+o 2 =so 2 +ho 2 

1.0E+12 

0.0 

503 


35 

hoso 2 +m=so 3 +h+m 

3.2E+16 

-0.81 

27024 



Ar,0.75/ H 2 ,2.5/ H 2 O,12.0/ 




e 

36 

hoso 2 +h=so 2 +h 2 o 

1.0E+12 

0.0 

0 


37 

hoso 2 +o=so 3 +oh 

5.0E+12 

0.0 

0 


38 

hoso 2 +oh=so 3 +h 2 o 

1.0E+12 

0.0 

0 


39 

hoso 2 +o 2 =so 3 +ho 2 

7.8E+11 

0.0 

330 



° Rate coefficients are in the form, k = AT n exp(- 0 /T). Units are K, cm, mol, and s. 

b k oo 

c ko 

a Falloff correction parameters. 
e Collision efficiencies. 

f 

Duplicated reaction. 
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Log(A i/cm mol s’ ) Log (k x /cm 3 mol 


H + 0 2 = OH + O 



10000 KIT 

Figure A 1(a). — Rate coefficients of Rl, H+0 2 — > OH+O. 



10000 K IT 

Figure Al(b). — Exploded plot of the Figure Al(a) for the temperature range of 950 to 1250 K. 
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Log (k 2 /cm 3 mol ' 1 s' 1 ) Log(A 2 /cm 3 mol' 


O + H 2 = OH + O 



10000 KIT 


Figure A2(a). — Rate coefficients of R2, 0+H 2 — > OH+H. 


o+h 2 = oh + o 



Figure A2(b). — Exploded plot of the Figure A2(a) for the temperature range of 830 to 1250 K. 
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Log(£ 21 ,o/cm mol 2 s 1 ) (k 9 y[M]/cm 6 mol 


h + o 2 + m = ho 2 + m 



o + so 2 +m= so 3 +m 

16.0 


15.5 


15.0 


14.5 


14.0 

0.0 10.0 20.0 30.0 40.0 50.0 

10000 KIT 

Figure A4(a). — Low pressure limit rate coefficients (k 2 1 , 0 ) of R21, CM-S0 2 +M — > S0 3 +M. 
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o + so 2 =so 3 



10000 K/T 


Figure A4(b). — High pressure limit rate coefficients (k 2] oo) of R2 1 , 0+S0 2 — > S0 3 . 
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Appendix B 


Hydroxysulfonyl (HOSO 2 ) 


A f fl°(298.15 K)— 391.2 ki mol 1 

Vibrational Frequencies (cm -1 ) and Degeneracies: 

3539.8(1) 

1309.2(1) 

1296.2(1) 

1097.2(1) 

759.3(1) 

544 (1) 

491 (1) 

437 (1) 

252 ( 1 )- torsion 

Ground State Quantum Weight: 2 ( 2 A’) 

Point Group: Ci 

External Rotational Symmetry Number: c=l 

Bond Length (angstrom): 

r(S-01)=1.612 

r(S-02)= 1.444 

r(S-03)=1.436 

r(Ol-H)=0.967 

Bond Angle (degree): 

0(S-O1-H)= 107.2 
0(O1-S-O2)=1O8.4 
0(O1-S-O3)=1O6.O 
t(02-S-01-H)=26.1 
t (03-S-0 1 — H)= 161.1 

Product of the Moment of Inertia: I x I y I z = 1.473 x 10 _1 14 g 3 cm 6 
Reduced Moment of Inertia: I r = 1.405x 10* 40 g cm 2 
Internal Rotational Symmetry Number: n=2 
Barrier to Internal Rotation: Vo=9.5 ki mol -1 

Enthalpy of Formation 

From Benson's estimate for D°298.i5(HO-S02)=36±3 kc al mol -1 ,l Af if (298.15 K) = — 409 ki 
mol -1 (-98 kca\ mol -1 ) was derived. This value of Af ff (298.15 K) gives A T H° (298.15 K )=6 
kca\ mol -1 for the reaction, HOSO 2 +O 2 — > HO 2 + SO 3 . 

Margitan's rate coefficients measurement for the above reaction put an upper limit for the 
reaction endothermicity, A r ff (298.15 K)<2 kca\ mol -1 (Z7°298.i5(HO-S02)<32 kc al mol -1 ). 2 
Gleason and Howard set a limit on the endothermicity of the reaction (1), A r H° (298.15 K)<3 
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Acal mol \ 3 which places a lower limit of Af// 5 (298.15 K)>-395.4 kJ mol 1 (>-94.5 Aeal mol 
'). Here we adopted Margitan's value of A r //° (298.15 K) = 2 kca\ mol -1 , which gives Af if 
(298.15 K)=-391. 2 kS mof 1 (-93.5 /teal mof 1 ). 

Heat Capacity and Entropy 

The molecular structure of HOSO 2 was taken from the results of MP2 level calculation with cc- 
pVQZ basis set by Majumdar et al. 4 The reduced moment of inertia, I r , for the hindered internal 
rotation was calculated by assuming HOSO 2 is essentially a cis form. 4 This assumption sets the 
internal rotational symmetry number, n=2. The rotational barrier was calculated using the 
equation, 2 too) /n=(Vo/2I r ) 1/2, where c is the speed of light and (ois the torsional frequency, 252 
cm -1 . 

Hashimoto et al. measured four vibrational frequencies of HOSO 2 in an argon matrix at 1 1 K. 5 
Nagase et al. confirmed the four vibrational frequencies using the same method of Hashimoto et 
al. and also calculated the geometry and vibrational frequencies at the unrestricted HF level with 
the 3-2 1G(*) basis set. 6 Later, Kuo et al. measured vibrational frequencies of HOSO 2 in the 
argon matrix at 12 K. 7 The four frequencies reported by Hashimoto et al. and Nagase et al. were 
again confirmed. An additional frequency corresponding to HOS bend at 1296.2 cm -1 was 
measured. 7 Here, we took five experimental frequencies (3539.8, 1309.2, 1296.2, 1097.2, 759.3) 
of Hashimoto et al., Nagase et al. and Kuo et al. and four computed frequencies (544, 491, 437, 
252) of Nagase et al. for the derivation of heat capacity and entropy. 
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Hydroxysulfinyl (HOSO) 


A f //°(0 K)— 237.4 ± 10 £J mol -1 

Vibrational Frequencies (cm -1 ) and Degeneracies: 

3476 (1) 

1351 (1) 

1073 (1) 

764 (1) 

343 (1) 

200 (1)- torsion 

Ground State Quantum Weight: 2 ( 2 A”) 

Point Group: C s 

External Rotational Symmetry Number: c=l 

Bond Length (angstrom): 

r(S-01)=1.661 

r(S-02)=1.482 

r(Ol-H)=0.983 

Bond Angle (degree): 

0(S-O1-H)=1O6.9 
0 (01-S-02)= 109.8 

Product of the Moment of Inertia: I x I y I z =2.539x 10 -115 g 3 cm 6 
Reduced Moment of Inertia: I r =1.585x 10" 40 g cm 2 
Internal Rotational Symmetry Number: n=2 
Barrier to Internal Rotation: Vo=7.9 kJ mol -1 

Enthalpy of Formation, Heat Capacity and Entropy 

The enthalpy of formation of HOSO and its geometry and vibrational frequencies have not been 
measured experimentally. The most reliable information is from theoretical calculations of 
Laakso et al. 1 Laakso et al. optimized geometry and vibrational frequencies at the 
MP2=FULL/6-31G* level. Energies were obtained with G2 theory. For the enthalpy of 
formation, the changes in G2 energy for the reaction, HOSO — » OH+SO (A T H°(0 K)=280.8 kJ 
mol -1 ) was obtained and then AfZ/°(0 K) of HOSO was derived using the measured Af7Z°(0 K) 
values of OH and SO in the JANAF tables. 2 Laakso et al. proved the reliability of their 
calculation by comparing their G2 level A r // , (0 K) for the reaction, HS+O 2 -+ OH+SO with that 
of the experimental value from the JANAF tables. Only a small difference in the reaction 
enthalpy, A r //° (0 K)= 0.3 kJ mol -1 resulted. However, a recent measurement and quantum 
calculation of D°(OH, 0 K) 3 ’ 4 gave A f //° (OH, 0 K)=8.85 ± 0.07 *cal mol -1 (37.03 ± 0.29 kJ mol" 
*). Thus the calculated A r Z/°(0 K) value for HS+O 2 — » OH+SO reaction by Laakso et al. would 
be different from 0.3 kJ mol -1 . Since D°(0 K) for HOSO dissociation from Laakso et al. is the 
only available value, we adopted this value and calculated the enthalpy of formation of HOSO 
using the newly measured Af/7° of OH and NASA value of SO at 0 K. The recommended 
enthalpy of formation of HOSO is (HOSO, 0 K)=A f B°(OH, 0 K)+A f /7°(SO, 0 K)- 
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279.2=37.03+4.7143-279.2 kJ mol" 1 ^237.4 + 10 kJ mol" 1 (-56.7 + 2.4 kca\ mol" 1 ) with an error 
bound given by Laakso et al. 

For the derivation of heat capacity and entropy we took the molecular geometry and vibrational 
frequencies of Laakso et al. The barrier to internal rotation (Vo=7.9 kJ mol" 1 ) was directly taken 
from Laakso et al.'s computed G2 barrier for cis to trans transition of HOSO. 
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Sulfinyl (HS0 2 ) 


A f /J° (298. 1 5 K)^l 78.2 ± 8.4 kJ mol" 1 

Vibrational Frequencies (cm" 1 ) and Degeneracies: 

2269(1) 

1600(1) 

1128(1) 

1065(1) 

810(1) 

431(1) 

Ground State Quantum Weight: 2 ( 2 A’) 

Point Group: C s 

External Rotational Symmetry Number: g=1 

Bond Length (angstrom): 

r(S-01)= 1.447 

r(S-02)= 1.447 

r(H-02)=1.381 

Bond Angle (degree): 

0(O1-S-O2)=123.2 
0 (H-S-02)= 105.6 
X (H-S-02(0 1 ))= 1 2 1 . 1 

Product of the Moment of Inertia: I x I y I z = 1 .485 x 1 0" 1 15 g 3 cm 6 
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Enthalpy of Formation 

Benson's estimate for the enthalpy of formation of HSO 2 was Af//°(298.15 K)— 177.4+ 16.7 ArJ 
mol -1 (-42.4 + 4 &cal moF 1 ). The corresponding £>°298 i5(H-S02)=41 Areal mol -1 was compatible 
with the fast rate coefficients measured in flame studies at T>1600 K for the three-body 
recombination reaction (H+S 02 +M=HS 02 +M).' Boyd et al. performed ab initio calculations for 
HS0 2 at the HF/STO-3G* level . 2 Using the calculated D 0 0 (H-SO 2 )= 142.1 ki moF 1 and the 
A fF/°(0 K) values of H and S0 2 in the JANAF tables , 3 A f /r(0 K)^220.4 kS moP 1 (-52.7 Areal 
mol -1 ) was derived for HSO 2 . Laakso et al. optimized geometry and vibrational frequencies for 
HSO 2 at the MP2=FULL/6-3 1G* level . 4 Energies were obtained with G2 theory. For the 
en thalpy of formation, the change in G 2 energy for the reaction, HSO 2 — » H+SO 2 ( A T H°(0 
K)=56.8± 10 ArJ mol -1 ) was obtained and then AfF/°(0 K) of HSO 2 was derived using the 
measured A f H°(0 K) values of H and S0 2 in the JANAF tables. A f //°(298.15 K) ofHS0 2 =- 
141.4 ± 10 ArJ mol _ 1 (-33.8 + 2.4 Areal moP 1 ) was recommended. Laakso et al. proved the 
reliability of their calculation by comparing their G2 level A r JJ°(0 K) for the reaction, HS+O 2 — > 
H+SO 2 with that of the experimental value from the JANAF tables. Considering the error limit in 
the: G2 A r /f°(0 K) value, the difference in the reaction enthalpy, A r F/°(0 K)=9.4 ArJ mol -1 was 
considered excellent. Recently, Denis and Ventura investigated the enthalpy of formation of 
XSO 2 , X=H, CH 3 using the B3LYP and B3PW91 methods with very large basis sets up to cc- 
pV 6 Z . 5 The enthalpy of formation of HSO 2 , AfF/°(298.15 K)— 178.2 ± 8.4 kJ mol -1 (-42.6 + 2 
foal moF 1 ) resulting from B3PW91 method was recommended. We adopted the enthalpy of 
formation of Denis and Ventura, Af7/°(298.15 K)— 178.2 ± 8.4 ArJ mol -1 for HSO 2 . 

Heat Capacity and Entropy 

Vibrational frequencies of HSO 2 have not been measured experimentally and only one 
theoretical calculation is available. For the derivation of heat capacity and entropy we used the 
vibrational frequencies calculated at the MP2=FULL/6-31G* level by Laakso et al . 4 The 
molecular geometry was taken from the theoretical calculation using the B3PW91 method with 
cc-pV 6 Z basis set by Denis and Ventura . 5 
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Mercapto (SH) 


A f fl°(0 K)= 1 33.03 JfcJ moF 1 
Enthalpy of Formation 

The current JANAF recommendations for Af/Z^O K)=136.49±5.0 ki moF 1 and Af// ) (298.15 
K)=l 39.33 ± 5.0 ki mol -1 are based on the measurements of HS ionization potential and the 
appearance potential of HS + from photoionization of H 2 S . 1 The previous NASA recommendation 
for A f //°(298.15 K) 2 was also 139.33 ki mol '. Continetti and Lee 3 performed photodissociation 
of H 2 S at 193.3 nm using H-atom photofragment-translational spectroscopy with mass 
spectroscopic detection. They found D°o(H-S)=349.27±2.89 ki moF 1 (83.48+0.69 A:cal moF 1 ) 
in the secondary photodissociation of HS radicals, formed in the primary H 2 S photodissociation 
experiments. Nicovitch et al. 4 measured the forward and reverse reaction rate coefficients of 
Br( 2 P 3 / 2 )+H 2 S=HS+HBr over the temperature range of 3 19 to 43 1 K. Using both the "second law 
method" and the "third law method," an average value of A r H° at a given temperature was 
obtained, then the values of A T H(0 K) and A r //°(298.15 K) were calculated by applying heat 
capacity corrections from data in JANAF tables for Br, H 2 S, HS and HBr. The reported values 
for the enthalpy of formation and the bond dissociation energy are Af/J°(0 K)=142.55±3.01 ki 
moF 1 (34.07 ±0.72 Acal moF 1 ), A f // > (298.15 K)=143.01 ±2.85 ki moF 1 (34.18±0.68 /fecal moF 
l ) and Z>° 0 (H-S)= r 348.23 ±3.26 ki moF 1 (83.23±0.78 /fecal moF 1 ), D°298.15(H-S)=351.96±3.10 
ki moF 1 (84.12±0.74 kcal moF 1 ), respectively. We adopted the bond dissociation energy of 
Continetti and Lee, £>°o(H-S)=349.27±2.89 ki moF 1 (83.48 ±0.69 Aval moF 1 ) and derived 
Af/J°(0 K)=Af//°o(S)+Af// ) o(H>-349.27=270.5126+21 1.8014-349.27=133.03 /fcJ moF 1 (31.80 
Acal moF 1 ). 
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Hydrogen Sulfur Oxide (HSO) 


Af//°(0 K>=-3.8 /fcJ moF 1 

Vibrational Frequencies (cm 1 ) and Degeneracies: 
2480(1) 

1164(1) 

1026(1) 
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Ground State Quantum Weight: 2 (2A”) 

Point Group: Cl 

External Rotational Symmetry Number: s=l 

Bond Length (angstrom): 

r(S-0)= 1.494 

r(S-H)=1.389 

Bond Angle (degree): 

0(H-S-O)= 106.6 

Product of the Moment of Inertia: I x I y I z =5.028x 10 -117 g 3 cm 6 


En thalpy of Formation 

Schurath et al. 1 measured the chemiluminescence spectrum of the X 2 A” <— A 2 A’ transition of 
HSO and DSO. The A r fl°(298.15 K) for the reaction HS+0 3 -» HS0+0 2 (HS+0 3 -> 

HS0*+0 2 , HSO* — > HSO) was estimated using the observed highest level of HS-0 vibration, 
v ' 3 =7, which corresponds to 19149 cm -1 (54.74 kcal mol -1 ). Using-A r H°(298.15 K)<54.74 
A:cal mol -1 and the values of the enthalpy of formation of HS and 0 3 available, Schurath et al. 1 set 
an upper limit of A ^(298. 15 K)< 14.9 Acal mol -1 . White and Gardiner 2 reanalyzed the results 
of Schurath et al. and presented 12.2 Af/7°(298.15 K)<14.1 kca\ mol -1 . Benson's estimate 3 was 
Aj77°(298.15 K)— 5 ± 4 Acal mol -1 Slagle et al. 4 suggested AfH° (298.15 K)— 3 kca\ mol -1 , 
based upon their kinetic work on the reaction 0( 3 P)+H 2 S — > HSO+H. Davidson et al. 5 studied 
the reactive scattering of O atoms with H 2 S molecules using a crossed molecular beam mass 
spectrometric detection of HSO radicals. The value of the enthalpy of formation obtained was 
Aj77°(298.15 K)— 1.4± 1.9 kca\ mol -1 (-6± 8 LJ mol -1 ). Balucani et al. 6 determined AfT/°(0 K) 
— 0.9 ± 0.7 kca\ mol -1 (-3.8 ± 2.9 kJ moT 1 ) from the analysis of high resolution crossed beam 
reactive scattering experiments on the reaction 0( 3 P)+H 2 S — > HSO+H. Here we took Balucani 
et al.'s value A f fl°(0 K)=-0.9± 0.7 *cal mol -1 (-3.8 + 2.9 kJ moT 1 ). 

Heat Capacity and Entropy 

The molecular geometries and vibrational frequencies of HSO and DSO radicals were 
determined by Schurath et al., 1 Kakimoto et al., 7 and Ohashi et al. 8 Following the measurements 
of Schurath et al., Kakimoto et al. and Ohashi et al. observed the doppler-limited dye laser 
excitation spectra of the A 2 A’(003)<— X 2 A”(000) vibronic transition of HSO and DSO. Here we 
took the structural parameters of HSO from those of Ohashi et al. and Kakimoto et al., instead of 
the "best fit" values used in vibronic band contour synthesis by Schurath et al. The product of 
moment of inertia, I x I y I z , was calculated using the directly measured rotational constants, B x , B y , 
and B z of Ohashi et al. 

The three ground state ( 2 A”) vibrational frequencies, H-SO stretching (coi), H-S-0 bending 
( to 2 ), and HS-0 stretching ( co 3 ), were presented by Schurath et al. 1 From v ’ 3 progression, w 3 
(1013 cm -1 ) was measured. The quantity co 2 was estimated (1063 cm -1 ) using the measured 
bending frequency of DSO and the isotope factor. For coi, the recommended — S-H group 
frequency (2570 cm -1 ) 9 was assigned. However, the harmonic force field analysis by Ohashi et 
al. 8 yielded 0) 2 =1 164 cm -1 and 0)3=1026 cm -1 . The force constant, k\, was estimated by applying 
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Badger's rule using the bond lengths and force constants of HS( FI ) and HS( X ) as references. 
From this estimated force constant o>i=227 1 cm -1 was calculated. 

Here we adopted w 2=1 164 cm" 1 and 0)3=1026 cm" 1 of Ohashi et al. We believe the 0)i value 
from — S-H group frequency is too high and that of Ohashi et al. is too low. Therefore we 
estimated the 0>i value in the following manner: Using the known H-O and H-OO stretching 
frequencies the force constant change from H-0 to H-OO was calculated. For H-OO stretching 
the force constant was reduced by 15% from that of H-O. Since HS and HSO are isovalent 
radicals of HO and HOO, respectively, we calculated the force constant of H-SO by applying the 
same 15% reduction to the H-S force constant. The coi value was calculated to be 
Wi=(1/27Ic)(&i/wh) 1/2= 2480 cm -1 , where c is the speed of light and wh is the mass of the 
terminal H atom. 
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